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We perform a numerical study of ghost condensation — in the so-called Overhauser channel — 
for SU(2) lattice gauge theory in minimal Landau gauge. The off-diagonal components of the 
momentum-space ghost propagator G cd (p) are evaluated for lattice volumes V = 8 4 , 12 4 , 16 4 , 
20 4 , 24 4 and for three values of the lattice coupling: (3 = 2.2, 2.3, 2.4. Our data show that the 
quantity <j> (p) = e bcd G cd (p) /2 is zero within error bars, being characterized by very large statistical 
fluctuations. On the contrary, |<^> 6 (p)| has relatively small error bars and behaves at small momenta 
as L~ 2 p~ z , where L is the lattice side in physical units and z fa 4. We argue that the large 
fluctuations for (f> b (p) come from spontaneous breaking of a global symmetry and are associated 
with ghost condensation. It may thus be necessary (in numerical simulations at finite volume) to 
consider |</> ! '(p)| instead of </> b (p), to avoid a null average due to tunneling between different broken 
vacua. Also, we show that </> 6 (p) is proportional to the Fourier-transformed gluon field components 
ATj(g). This explains the L~ 2 dependence of \<f> (p)[, as induced by the behavior of \A b ^(q)\- We 
fit our data for j</> b (p)| to the theoretical prediction (r/L 2 + v)/(p 4 + v 2 ), obtaining for the ghost 
condensate v an upper bound of about 0.058 GeV 2 . In order to check if v is nonzero in the continuum 
limit, one probably needs numerical simulations at much larger physical volumes than the ones we 
consider. As a by-product of our analysis, we perform a careful study of the color structure of the 
inverse Faddeev-Popov matrix in momentum space. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.Gc, 14.80.-j 



INTRODUCTION 



The QCD vacuum is known to be highly non-trivial 
at low energies [lj. This non-trivial structure manifests 
itself through the appearance of vacuum condensates, 
i.e. vacuum expectation values of certain local operators. 
In perturbation theory these condensates vanish, but in 
the SVZ-sum-rule approach |2|, |3j they are included as a 
parametrization of non-perturbative effects in the eval- 
uation of phenomenological quantities. Of course, the 
local operators considered must be gauge invariant in or- 
der to be added as higher-dimension operators (higher 
twists) to the usual perturbative expansion of physical 
quantities. The two main such operators are a s F tlv F t11 ' 
and m q ip q ij) q ; their vacuum expectation values are the 
so-called gluon and quark condensates. Both of these 
operators have mass dimension four. 

In recent years, (gauge-dependent) condensates of 
mass dimension two have also received considerable at- 
tention Hj EBJl jlMMEE B1I|JI1JI11I| 



Re- 



a physical interpretation to the condensate (A^A 1 ^) 
cently, evidence that the expectation value ( A b A b ) may 
be gauge independent was presented in 

m ee m est 

At the same time, the authors of Refs. |4J,|5| have shown 
that the minimal value of (A b A b ) (on each gauge orbit) 
has a gauge-invariant meaning. Let us note that this 
minimal value is achieved by considering absolute min- 
ima in the minimal Landau gauge, i.e. for configurations 
belonging to the so-called (Landau) fundamental modu- 
lar region [38| . 

Other vacuum condensates of mass dimension two 
considered by several groups are the ghost condensates 




In particular, the gauge condensate (A b A b ) has been 
largely studied, since it should be sensitive to topological 
structures such as monopoles 0,0- Thus, this parameter 
could play an important role in the quark-confinement 
scenario through monopole condensation |35l 1-3(1 l37| . 
Moreover, the existence of a gauge condensate would im- 
ply a dynamical mass generation for the gluon and ghost 
fields. This result has already been obtained in various 
gauges [S 113 El |M IH, but clearly, since the opera- 
tor A b A b is not gauge invariant, it is difficult to assign 



(f bcd rj c r] d ), (f bcd r) c r] d ) and (f bcd r] c rj d ) . These conden- 
sates were first introduced in SU( 2) gaug e theory in max- 
imally Abelian gauge (MAG) E El El iM III El 
More recently, the same condensates have been studied 
in other gauges [1 E U El El E3 , such as the Curci- 
Ferrari and the Landau gauges. In all cases it was found 
that the ghost condensates are related to the breakdown 
of a global SL(2, R) symmetry [rlE^. In MAG the diag- 
onal and off-diagonal components of the ghost propaga- 
tors arc modified 0, E3 by ghost condensation. Similar 
results have been obtained in other gauges 13 y(j • In 
particular, in Landau gauge it was shown [13j that the 
off-diagonal (anti-symmetric) components of the ghost 
propagator G cd (p) — or, equivalently, the quantity <j> b (p) 
- are proportional to the ghost condensate v. In the 
same reference it is argued that the breaking of the global 
SU(N C ) color symmetry, related to this non-zero expec- 
tation value for G cd (p), may occur only in the unphysical 
sector of the Hilbert space and therefore should not be 
(explicitly) observable in the physical subspace. Very re- 



2 



cently [23|, the effects of the ghost condensate (f bcd r] c r] d ) 
and of the gauge condensate (A b ^A b ^) have been consid- 
ered together for the SU (2) case in Landau gauge, show- 
ing that ghost condensation induces a splitting in the 
gluon mass related to the gauge condensate. 

It is interesting that the existence of different ghost 
condensates has an analogue in the theory of supercon- 
ductivity 0. In particular, the condensates (f bcd r] c r] d ) 
and (f bcd T] c r] d ) correspond to the so-called BCS effect 
|43l 23 j m which particle-particle and hole-hole pairing 
occur, while the condensate (f bcd r] c ri d } is analogous to 
the Overhauser effect (i^, in which particle-hole pair- 
ing occurs. Also, in the Curci-Ferrari gauge, both effects 
have to be considered |13( in order to obtain an action 
invariant under the so-called Nakanishi-Ojima (NO) al- 
gebra. It is this NO algebra that includes an SX(2, R) 
sub-algebra. 

Finally, a mixed gluon-ghost condensate of mass di- 
mension two has also been studied by several authors 
PUEEHHUEJ], using various gauges. This mixed 
condensate is of particular interest when considering in- 
terpolating gauges 0|- Indeed, it allows to generalize 
and relate results obtained in different gauges for the 
gauge condensate (A'^A'p and the ghost condensates. 
Moreover, in MAG this mixed condensate would induce a 
dynamic mass for the off-diagonal gluons |20| , giving sup- 
port to the Abelian-dominance scenario pa l47t l48l |49| . 
Thus, the various gauge and ghost condensates could all 
play an important role in the dual superconducting sce- 
nario of quark confinement |50l ISlL 152 ] , being related to 
monopole condensation and to Abelian dominance. 

Possible effects of the gauge condensate (A b ^ A b ) on 
propagators and vertices (in Landau gauge) have been 
studied thr ough lattice simulations in Refs. [5^, E^, 
IH |H |H H3, yielding {A^AD k 3 GeV 2 . On the 
other hand, the ghost condensate has been considered 
numerically only in Ref. |60j . for minimal Landau gauge 
in the Overhauser channel, both for the SU (2) and SU (3) 
cases. The authors find that the color off-diagonal anti- 
symmetric components of the ghost propagators are iden- 
tically zero, while the symmetric components are zero 
within error bars, but have fluctuations proportional to 
A~ 2 p~ 4 , where p is the lattice momentum and N is the 
lattice side. 

In this paper we carry out a thorough investigation of 
ghost condensation in the Overhauser channel for pure 
SU(2) Yang-Mills theory in minimal Landau gauge. In 
particular, we evaluate numerically the off-diagonal com- 
ponents of the ghost propagator G cd {p) as a function of 
the momentum p. We find that (<j) b (p)) = e bcd (G cd {p))/2 
is zero within error bars, but with large fluctuations. On 
the other hand, if we consider the absolute value of 4> b (p) , 
we find that ( | <j) b (p) \ ) behaves at small momenta as 
L~ 2 p~ z , where L and p are, respectively, the lattice side 
and the momentum in physical units and z ~ 4. As we 
will show [see Eq. I|t)5|) ]. the sign of 4> b (p) is related to the 



sign of the Fourier-transformed gluon field components 
A(q). At the same time, one should recall that a nonzero 
value for 4> b {p) is related to the spontaneous breaking 
of a global SL(2,R) symmetry Thus > the situa " 

tion here is similar to the one encountered in numerical 
studies of spin systems with nonzero spontaneous mag- 
netization, in the absence of an external magnetic field. 
Indeed, since spontaneous symmetry breaking can occur 
only in the thermodynamic limit, when performing nu- 
merical simulations at finite volume one always finds (af- 
ter sufficiently many Monte Carlo st eps) a null magneti- 
zation for all nonzero temperatures [102| . This is due to 
the fact that, at finite volume, the probability that the 
system may pass from a state near a broken vacuum to a 
state near a different broken vacuum is always nonzero. 
To avoid this problem, the usual solution is to consider 
the absolute value of the magnetization or the root mean 
square order parameter, i.e. ( | M \ ) or ( M 2 ) 1 / 2 . The re- 
sults obtained using these two quantities do not agree at 
finite volume, but they produce the same infinite- volume- 
limit result |6l], Fig. 2.13]. Alternatively, one can also 
simulate spin models with an external magnetic field H 
[63| and then consider the limit H — > 0. In our case, 
one could set the sign of 4> b {p) by adding to the action a 
constant external chromo- magnetic field coupled to the 
gluon field or, equivalently, one can fix this sign by a 
global gauge transformation. This is analogous to using 
a global rotation in a spin system in order to fix a direc- 
tion for the magnetization. In all cases, i.e. introducing 
an external magnetic field, making a global transforma- 
tion, or taking the absolute value of the order parameter, 
one explicitly breaks the global symmetry of the system 
choosing one of the (equivalent) possible vacua. In the 
present work we will consider the absolute value of <j> b (p) , 
but we also describe (in Section Hll Aj l a possible global 
gauge transformation to select a unique broken vacuum. 

Another indication of spontaneous symmetry break- 
ing is the direct investigation of the statistical distribu- 
tion of the order parameter, obtained from histograms 
of the produced Monte Carlo data. Let us recall that, 
given a random variable x with a Gaussian distribu- 
tion e~ x I ( 2a ) / \/2~7Ta, corresponding to a null mean value 
and a standard deviation ^/a, one finds that the vari- 
able | a; | has mean value yj2a/ir and standard deviation 
^{-k —2) a/n w 0.75 (|x|). Thus, if one evaluates the 
quantity x and finds that its average is null within er- 
ror bars while \x\ has a relatively small error, one should 
conclude that the data for x do not correspond to a Gaus- 
sian distribution and that the fluctuations observed are 
not just statistical ones, since in this case they should 
be large also when considering the absolute value of x. 
This is indeed what happens for the magnetization of 
a spin system in the broken (low-temperature) phase, 
as can be easily verified by comparing a histogram of 
the data with a Gaussian distribution. This is true 
both for discrete symmetries — as in the Ising model, 
for which the distribution of the magnetization has two 
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peaks |6ll l62j ] — and for continuous symmetries, such 
as for the three-dimensional 0(iV)-symmetric nonlinear 
cr-models. In this case (see Fig. QJ, the shape of the sta- 
tistical distribution of the data is similar to a Gaussian, 
but with a larger width around the peak of the curve. 
Moreover, the width becomes larger as the temperature 
of the system decreases. On the contrary, in the unbro- 
ken phase, the histogram of the data is indeed a Gaussian 
distribution (see again Fig.QJ. Note that a simple way of 
evaluating the Gaussian character of a distribution 
is to consider the so-called Binder cumulant 



for the order parameter x. In the infinite-volume limit 
this quantity is zero if x has a Gaussian distribution, i.e. 
in the unbroken phase, and is equal to 2/3 if the distri- 
bution of x has two peaks, as in the broken phase of the 
Ising model. Thus, in general one should expect U 7^ 
in the broken phase, where the value of U depends on 
the distribution of x. In the case of the ghost condensate, 
based on the discussion above, we argue that the occur- 
rence of large fluctuations for <fr (p) but not for | <fi b (p) | 
implies that the global symmetry SL(2, R) is indeed bro- 
ken. We will analyze the statistical distribution for <fi b (p) 
111 Section |V| 

As mentioned above, we show in Eq. (|65f) that the 
quantity <fi b (p) is proportional to the Fourier-transformed 



gluon field A(q). This dependence can explain the 1/L 2 
behavior of (j) b (p) as induced by A(q). Indeed, for q = it 
was proven 64] that ( | A(Q) | ) goes to zero in the infinite- 
volume limit at least as fast as 1/L. Furthermore, in 
Section [V] we verify from our data that | A b (q) | behaves 
as 1/L 2 . We stress that this result does not imply that 
|0 b (p)| is zero at infinite volume, because (\4> b (p)\) is 
not simply proportional to ( | A b (q) | ), as can be seen in 
Eq. 

The paper is organized as follows. In Section [H] we 
review the definition of the Faddeev-Popov operator — 
from which we obtain the ghost propagator — on the 
lattice (for the minimal Landau gauge) and in the con- 
tinuum (in Landau gauge) for a general SU(N C ) gauge 
group. In Section ILTT1 we define <fi b (p) and show that it is 
proportional to the Fourier-transformed gluon field A(q). 
In that Section we also discuss how to find a global gauge 
transformation that changes the sign of cf> b (p). Our nu- 
merical simulations are explained in Section ITvl and the 
results are reported in Section [V] Let us stress that, 
as a by-product of our analysis, we also present a careful 
study of the color structure of the inverse Faddeev-Popov 
matrix in momentum space. At the same time, we con- 
sider the infinite-volume limit for the momentum-space 
gluon field components A(q). Finally, in Section IVII we 
draw our conclusions. 



II. THE LANDAU FADDEEV-POPOV 
OPERATOR 
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FIG. 1: Histograms for the first component of the magneti- 
zation in the 3d 0(4) nonlinear er-model. Solid, dashed and 
dotted lines correspond respectively to inverse temperature 
/3 = 0.9 (in the symmetric phase), /3 — 0.98 and /3 = 1.2 
(both in the broken phase). In the three cases we consid- 
ered 100,000 independent configurations. For comparison, we 
show (dot-dashed line) a Gaussian with the same mean and 
variance as the case /3 = 0.98. 



Let us recall that, on a d-dimensional lattice, the (min- 
imal) Landau gauge is obtained, in the SU (N c ) case, by 
considering the functional 

£u[9] = "dVEE^f [g{x)U,{x)g\x + e,)\ (2) 

X fl ° 

= -^?EE JJf [9{x)Uv(x)g*(x + eJ 

X fl ° 

+ g{x + e ll )Ul(x)g\x)] . (3) 

Here {U^(x)} is a given (i.e. fixed) thermalized lattice 
configuration, V is the lattice volume and the gauge- 
fixing condition is imposed by finding a gauge transfor- 
mation {g(x)} that brings this functional to a minimum. 
Both the link variables U^(x) and the gauge transforma- 
tion matrices g(x) are elements of the SU(N C ) group (in 
the fundamental N c x N c representation). Also, 1L indi- 
cates the identity matrix and e M is a positive unit vector 
in the [i direction. Note that we set the lattice spacing a 
equal to 1 in order to simplify the notation. 

Let us also consider [2a E5 a one-parameter subgroup 
g(r;x) — exp [ i r r y b (x)X b ] , where r is a real param- 
eter, A b are the N 2 — 1 traceless Hermitian generators 
of SU(N C ) (also in the fundamental representation), the 
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components of 7 b (x) are real and the sum over the color 
index b is understood. We consider generators A b nor- 
malized such that Tr(A b A c ) = 2S bc . In the SU(2) case 
these generators are the three Pauli matrices a b . Using 
this one-parameter subgroup we can regard £jj as a func- 
tion of the parameter r. Its first derivative with respect 
to r calculated at r = is given by 



If {U^ (x)} is a stationary point of £(t) at r = then 
we have £ (0) = for every {^ b (x)}. This implies [see 
Eq. Q] 



(V • A b ) (x) = 



(15) 



£(0) 



for any x and &, which is the lattice formulation of the 
usual Landau gauge-fixing condition in the continuum. 



dVN c 



^2 l h ( x ) ( v ' Ah ) ( x ) ■ ( 4 ) Let us recall that the continuum gluon field A^\x) [T 



Here 



(y-A b ){x) =Y J Al{x)-Al{x-e,) (5) 



is the lattice divergence of the gluon field A b (x) , defined 

as 



4(z) = ±Tr[A fl (x)\ b ] 



with 110 



Mx) = i [U^x)-Vl(x) 



M ' Jtraceless 



(G) 



(7) 



or, equivalcntly, 



.!,,(.,•) = ^ [U^x)-Ul(xY 



•1L- 



Tr 

2i7V, 



^(x)-^)] . (8) 



The last term in the above equation is zero in the SU (2) 
case and the definition of the lattice gluon field simplifies 
to 



1 



A^x) = - [U^(x)-Ul(x)] . 
Then, one can write the SU{2) link variables U^{x) as 



(9) 



with 



U^x) = *7°(x)i_ + iA b {x)a b (10) 



[u°M 2 + E [4(*)] 2 = 1 ■ (") 



Also, note that, since the generators X b are traceless, the 
term proportional to the identity matrix 1L in Eq. (JSJ 
does not contribute to A b (x) [see Eq. ©]. Thus, for any 
N c > 2 we find 



.!„(,) = AM;;.,) = ^- Tr [Apix) A 6 ] (12) 



A^ 
A 6 



Tr{[U,(x)-Ul(x)] \ b } (13) 
ImTr [[/^(x) A 6 ] . (14) 



is related to the link variables U^x) by the relation 



exp 



i a g A<£>(x) 



(16) 



where #n is the lattice bare coupling constant. This 
yields jto^ 



A^x) = ag Af{x) + 0(a 3 ) 
Finally, using 

V M = aS M + 0(a 2 ) 
we obtain from Eq. I|15|) the result 

(d ■ A b ) (x) = 0{a) . 



(17) 



(18) 



(19) 



Still considering the one-parameter subgroup g(r; x) 
and the general SU (N c ) case we can evaluate the second 
derivative of £ (r) at r = 0. We obtain 



£ (0) = 



dVN, 



dVN L 



J2 E^^wMfM (20) 



x,y b,c 



£$> b (z) (M u7 ) b (x), (21) 



with |36 



(M u7 ) b (x) = X; i t( a; )[7 c (a:)-7 c (a: + ^)] 

A" 

+ F>-e M )[ 7 c (x) - 7 C (^-^)] 
+/ Mc [A^^ix + e^) 

-Afa-eJy'ix-eJ] . (22) 

Here A4[f(x,j/) is the lattice Faddeev-Popov matrix and 
f bcd are the (anti-symmetric) structure constants defined 
by [A b ,A c ] = 2if bcd X d . We also have 



Tr 



X b A c I U»(x) + UUx) 



2 ' 2 



(23) 



In the SU{2) case one finds r bc (x) = S bc Tr C/ M (x)/2 
and / 6od = e bcd , where e bcd is the completely anti- 
symmetric tensor. In this case, Eqs. Il'ill) and l|22|l agree 
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with Eq. (11) of Ref. Note that, by using Eqs. © 

and l|15fl . we can re- write Eq. I)22|) as 

(M u7 ) h (x) = ]Tr^(x)[ 7 c (x) - 7 ^ + e M )] 

+ r^-e M )[ 7 c (x) - 7 c (x-eJ] 

+ f bdc {A d (x) [ 7 c (x + e p ) + 7 c (x)] 

-Afa-eJl-fix-eJ + 7 c (x)]} .(24) 
We then obtain 
(Muj) b (x) = J2 ( D ^) b (*) " (^7) b (a: - e p ) (25) 

= (V • D^) b (x) , (26) 

where we used the definition (0) of the lattice divergence 
and 

{D^f{x) = r b ;(x)[ 7 c (x)~^(x + e,)} 

+ f bdc A d (x) [ 7 c (x + e M ) + 7 c (x)] (27) 

is the (negative of the) lattice covariant derivative |3S| 
applied to j c (x). Equivalently, we can write Eq. (|22|l as 

(Mujfix) = £lt(z)[ 7 c (z) -^(x + e^)] 

+ r>-e p )[ 7 c (x) - 7 c (x-e M )] 
+ / Mc {^( a: ) [ 7 c (x + e M )- 7 c (x)] 

- Afa-eJtfix-eJ - 7 c (x)]} (28) 

and it is evident that the second derivative £ (0) is null 
if the vector 7 c (x) is constant, i.e. the Faddeev-Popov 
matrix has a trivial null eigenvalue corresponding to a 
constant eigenvector. It is also clear that, if {U^ (x)} is a 
local minimum of £ (r) at r = 0, then the matrix Mu is 
positive definite (in the subspace orthogonal to the space 
of constant vectors). 

Let us recall that, in the continuum, if one considers 
a gauge-fixing condition [F(A 9 )] b (x) =0 and the gauge 
transformation g(x) — e lg ° e ( x ' x , the Faddeev-Popov 
operator can be determined using the expansion 

[F(M)] b (x) = [F(A)} b (x) 

+ f d 4 yJ2M bc (x iy )e c (y) + ... .(29) 

In the Landau case one has [F(A)] b (x) = (d^ A b J(x), 
giving 

M hc {x,y) = -d^D h °{x) 8{x-y) (30) 
= -D b ;(x)d,6(x-y) (31) 
= - [S bc - 2 50 / Mc A* (x)] d^6(x- y) (32) 
= - [S bc d„ + 2 g Q f bcd A d (x)] d^ix-y), (33) 



where D h ^{x) is the (continuum) covariant derivative. 
Note that, with our normalization for the generators X b , 
we have that 2 A® (x) is the standard continuum gluon 
field. 

Clearly, we can expand Eq. (|16[) in the limit a — > as 

C/°(x) = 1L + ioflb^^) + °(« 2 )' ( 34 ) 

yielding 

r^ c (x) « ^ c + 0(a 2 ). (35) 

Also, after Taylor expanding the vectors 7 & (x ± e M ) = 
7 b (x ± ae^) to order a 2 in Eq. (|22p. one obtains 

a^M^y) = M bc (x,y) + 0(a) , (36) 

with A^ bc (x, y) given in Eq. I|32|) . At the same time, from 
Eq. H27fl we find the relation between the lattice and the 
continuum expressions for the covariant derivative 

D b ;(x) = -aD^ b %x)+0(ai) (37) 

and it is immediate to obtain Eq. i|5U|) . using Eqs. i|2l)|) 
and HH}. 

III. THE GHOST CONDENSATE 

Let us note that the lattice Faddeev-Popov matrix 
Mjj{x,y) is symmetric with respect to the exc hang e 
(b, x) <-» (c,y). Indeed, from Eq. 1)221) we can write [7(j 

M b u c (x,y) = " Wv] 

A* 

+ (x — e^) [(5x,a ~ <^x— e M ,i/] 

-^(x- e/1 )«5 x _ e ^] (38) 

and it is straightforward to verify that A^[f(x, y) = 
A^^(y, x). Moreover, using the fact that the inverse of 
a symmetric matrix is also symmetric, we have that the 
same property is satisfied by the inverse matrix 

G bc (x,y) = (Mu x ) bc (x,y) . (39) 

Let us recall that, if a matrix M bc (x, y) has two different 
sets of indices, it can always be decomposed as 

M bc (x,y) = M bc s (x,y) + M b <(x,y) 

+ M br s (x,y) + M bc a (x,y), (40) 

where the sub-scripts s and a indicate symmetry and 
anti-symmetry with respect to the color or to the space- 
time indices (in this order). Furthermore, if the ma- 
trix M (x, y) is symmetric with respect to the exchange 
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(b,x) <-> (c,y), as in the case of the Landau Faddeev- 
Popov matrix Mjj, then this decomposition simplifies to 



M bc (x,y) = M bc s (x,y) + M b a c a (x,y) . 



(41) 



Indeed, for the Faddeev-Popov matrix M.\j(x,y) [see Eq. 
(|38[l ] we have that the terms containing Y b f{x) are sym- 
metric under the exchange b <-> c and under the exchange 
x <-» y. At the same time, the terms containing f bcd are 
anti-symmetric in the color indices and in the space-time 
coordinates x and y. Thus, in the SU(2) case, we can 
write [see Eqs. |22J ancl 

A4£(a:,y) = S bc S(x,y) ~ f bcd A d {x,y) (42) 

with S(x,y) — S(y,x) and „4 d (x,y) = ~A d (y,x), corre- 
sponding to the decomposition (|41|) . Note that the above 
formula also holds in the continuum [see Eq. for any 
N c > 2. Then, the inverse Faddeev-Popov matrix is given 
by 

(Mu l ) bc (x,y) = {5 be S(x,w)[S ec 6 w , v 

- f^S-^w^A^y)]}- 1 , (43) 
where the sum over repeated indices is understood. This 



implies that the inverse matrix G 



defined in Eq. 



(44) 
(45) 

(46) 



above can be evaluated using |106 

G bc (x,y) = S bc S-\x,y) 

+ f bcd S- 1 {x,z)A d {z,w)S- 1 {w,y) 
+ [f bdh f dc] S-\x,z)A h {z, S )S~\ S ,r) 
A'ir^S^iw^)] + ... 

We can also write 

G bc (x,y) = G bc (x,y) + G bc (x,y). 

Here 

G b f(x,y) = S^S-^y) + f bdh f dci [ 

S- 1 A h S- 1 A j S~ 1 ] (x,y) 

includes all terms of Eq. (|44|l with an even number of 
factors A d (z,w) and 

G b c (x,y) = f bcd [S- 1 A d S- 1 ] (x,y) + f d h f de jf ed { 

S- 1 A h S- 1 A j S- 1 A 1 S- 1 ] (x,y) + ... (47) 

is obtained by considering all terms with an odd number 
of factors A d (z,w). Then, it is easy to verify that 

1 

2 

G b e c (x,y) 

\ [G bc {x,y) ~ G cb (y,x)' 
G b c (x,y) , 



G b s c s (x,y) = - [G bc (x,y) + G cb (y,x) 



G b a c a (x,y) 



(48) 
(49) 
(50) 
(51) 



[whe re we considered again the decomposition Q41J1] . and 
that [l07| 



G b c (x,y) = f bde S-\x,z)A e (z,w)G dc (w,y) (52) 
G b e c (x,y) = S^S-'&y) 

+ f bdh S- 1 (x,z)A h (z,w)G d c (w,y). (53) 



Thus, for the ghost propagator G(x,y) — defined on the 



lattice as Tr Mj} /(N^ — 1) — we obtain 



G ( x ^ = W^l GbC{x>y) = T^G b e c (x,y), (54) 



N?-l 



while for ghost condensation (in the Overhauser channel) 
we consider |l3j the quantity 

<p\x,y) = ±f b ° d G cd (x,y) = ±f bcd G cd (x,y) (55) 



bed fceh 



f bcd f 



Also, in momentum space we have 

§ bc 



■S- 1 (x,z)A h (z,w)G e e d (w,y). (56) 



G(P) 



V(Ni-l) 



-Y,G b e c {x,y)< 



-ip{x-y) 



(57) 



x V 



V(N* 



— Y,G b e c (x,y)cos{p(x-y)}, (58) 



x y 



where in the last equation we used the symmetry of 
G bc (x,y) with respect to the space-time coordinates x 
and y. At the same time, we find 



jbed 

Vn~ c 
. / 



- ip (x—y) 



(59) 



xy 
bed 



J2Gf(x,y)sm[p(x-y)}. (60) 



■r V 



Considering the result (|56|l we can also write 



fbcd fceh 
\p) = 1 1 /-I o-l 



p\S-'\ q ){q\A h \k)(k\Gt d \p), 

(61) 

where we used Dirac notation and sums over q and k are 
understood. In order to simplify our analysis we can work 
in the continuum. Then, the matrix S(x,y) is given by 
-d 2 x 6(x - y) and A h = 2 g S(x - y) [see Eq. P5|) ]. 

Thus, we obtain 



and 



(plS-'lq) = 6(p-q)p- 2 (62) 



(q\A h \k) = 2igoA h (q-k)k^. (63) 
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Here A 1 ^ (p) is the Fourier-transformed gluon field defined 
(in the continuum and in d dimensions) as 



( 27r ) X 



x) e 



-i p X 



(64) 



Using the above results and ( k \ G ed \p) = G ed {k,p) we 
find 



fbcdjceh 2tgok f < A*(p-k)Gf(k,p) 
N c p 2 



(65) 



where a sum over the momenta k is understood. Note 
that, if G bc (k,p) is diagonal in color space, i.e. if 

G b e c (k lP ) = S bc G bb (k, P ), 



Eq. H65|) simplifies to 



Hp) 



f bce f ceh 2tg k^ A } l{p-k)G d e d {k,p) 



Nr. 



2ig k f , A^p-k) G d /(k,p) 



(66) 

(67) 
(68) 



where we used f bce f ceh = N c 6 bh . 

Finally, from Eq. (|53|> one obtains that the off-diagonal 
symmetric components of the ghost propagator G bc (x,y) 
are proportional to the off-diagonal anti- symmetric com- 
ponents G bc (x,y), being 

G b e c (x,y) = f hdh S- 1 {x,z)A h {z,w)G dc {w,y). (69) 

In analogy with Eq. i|65|) , one finds (in momentum space 
and in the continuum) 



G b e c ( P ) = I 



h,Ui 



2ig k fi A h (p~k)G d c (k,p) 



P 



(70) 



At the same time, for the ghost propagator, we can write 
[using Eqs. (JSU and J5D] 

G(x,y) = S^ix^) 



f 



cdh 



+ ^T^J S ^(^)A h (z,w)G d c (w,y) (71) 
S^ 1 (x,z) 5{z-y) 



A 



N, 



(72) 



Note that the above expression can be used to derive Eq. 
(27) of Ref. HU- 



A. Global gauge transformation 

Let us note that the term A b (x,y) is linear in the gluon 
field A d (x) [see Eqs. and gSJl]. Therefore, if we 
apply the transformation 



A d (x)^-A d (x), 



(73) 



then the ghost propagator G(x,y) does not change, while 
the quantity <f) b {x,y) gets multiplied by —1. Moreover, 
as we have seen in the previous section, when G bc (x,y) is 
diagonal in color space, the quantity 4> b {p) (in the con- 
tinuum) is proportional to the Fourier-transformed gluon 
field A^ [see Eq. ©]. 

Let us also stress that, for a given color index b, one can 
always find a global gauge transformation that changes 
the sign of A h {q). For example, in the SU(2) case 
and with 6 = 3 this can be achieved by using the x- 
independent gauge transformation 



g = i a 1 5i + ia 2 g 2 



(74) 



with g\ + g\ — 1. Indeed, if g = i J^b 0-6 9b> then the 
quantity 



1 

V' 



- iqx 



,t 



.9 



U°Jq)L + z<T c Al(q) 



becomes 

U°Jq)l. + 2iA c Jq)g c a b g b - iA c Jq)a c 



(75) 
(76) 

(77) 



Thus, if (?3 = one immediately finds that the sign of 
A^ (q) has changed (for any value of the index p and any 
momentum q). Here we considered the decomposition 
(|10fl for the SU{2) link variables U^(x) and the defini- 
tions 



% q x 



(78) 
(79) 



In the general SU(N C ) case (with N c > 2), changing 
the sign of A b (q) gives the relation 



Im Tr 



= — Im Tr 



£ U^x)e~ ldx \ b 



or 



Im Tr 



J2u,(x)e-^* (g^\ b g + X b ) 



(80) 



0, (81) 



where g is a global gauge transformation and we used 
Eq. I|14|) . Then, a sufficient condition to satisfy the above 
equation is given by 



gi\ b g + X b = 0. 



(82) 
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Lattice Volume V 


8 4 


f2 4 


f6 4 


20 4 


24 4 


No. of Configurations 


fOOO 


700 


400 


300 


too 



TABLE I: Lattice volumes and total number of configurations 
considered. 



p 


2.2 


2.3 


2.4 


a a 2 


0.220(9) 


0.136(2) 


0.071(1) 


a(GeV- r ) 


1.07(2) 


0.838(6) 


0.606(4) 



TABLE II: Lattice string tension and lattice spacing (in phys- 
ical units) for each lattice coupling /3 considered in our simu- 
lations. We used yfa = 0.44 GeV as physical input. 



implying 



{\\g} = 0. 



(83) 



Clearly, in the SU(2) case with 6 = 3 this yields g 
ia b gb and g% = 0. Note that the choice 



<Ti4(<?) + <T 2 Al{q) 



(84) 



Am 



does not to change the values of A^q) and A 2 (q). Of 
course, any other choice for g\ and g 2 (with g^ — 0) would 
still change the sign of (q) but would modify the values 

of the other two color components of A^q). In any case, 
a global gauge transformation does not change the value 
of U®(q) in Eq. (|76J) above, implying that the quantity 

r ~ i 2 

S c A c Aq) is also left unchanged. 



IV. LATTICE SIMULATIONS 

We performed numerical simulations of pure SU (2) 
Yang-Mills theory in Landau gauge consi derin g the stan- 
dard Wilson action in four dimensions 108J. To ther- 
malize the field configurations we used the heat-bath al- 
gorithm accelerated by hybrid overrelaxation (for details 
see jzHz3lz3)- In Tabled we show the lattice volumes 
used in the simulations and the corresponding numbers 
of configurations. For each volume we have considered 
three values of the lattice coupling: f3 = 2.2, 2.3 and 2.4. 
The corresponding string tensions in lattice units [T^, 13 
are given in Table ILT1 In the same table we also report 
the lattice spacing in physical units for each coupling j3, 
using y/a = 0.44 GeV as physical input. 

The minimal (lattice) Landau gauge is implemented 
using the stochastic overrelaxation algorithm [zS Eil Iz3 ■ 
We stop the gauge fixing when the quantity 



(V-AY 



(W-A b ) (*)]- 



(85) 



is smaller than 10 -13 . In order to check for possible 
Gribov-copy effects we have also done (for each ther- 
malized configuration) a second gauge fixin g u sing the 
so-called smearing method [78|. (See again |73| for de- 
tails.) 

Ghost condensation in the Overhauser channel is stud- 
ied by evaluating (here N c — 2) 



(i4> b (p)) 



^bcd 

TV 



J2(G cd (x,y)) S m[p(x-y)} (86) 



as a function of the lattice momentum p, with G cd (x,y) 
defined in Eq. l|39|l . For the inversion of the lattice 
Faddeev-Popov matrix Mjj we employ a conjugate gradi- 
ent (CG) method with even/odd preconditioning. Since 
this inversion has to be performed in the sub-space or- 
thogonal to the constant modes — corresponding to null 
eigenvalues of M.\j — we verify at each CG iteration if 
the solution has indeed a zero constant component. 

We also consider the average of the absolute value of 
4> b (p). Using the continuum result Q68[l. we should obtain 



'(P) 



650 kfj 



Ai(p- 



k) 



Gf(k,p)] 



P 



(87) 



This result is valid when the propagator G e e d {k,p) is diag- 
onal in color space, but it should give a good approxima- 
tion for ( J 4> b (p) J ) also if Gl d (k,p) is strongly diagonally- 
dominant. Let us recall that a mat rix My is strongly (or 
strictly) diagonally-dominant |l09j ] if Ma > YljjU Mjj for 
any i. We will verify this hypothesis in Section fVl 

Notice that a generic lattice momentum p has com- 
ponents (in lattice units) p^ — 2 sin (ttp^ a/L^) and 
that the magnitude (squared) of the lattice momentum 
is given by p 2 = p 2 . Here = aN^ is the physical 
size of the lattice in the [i direction, the quantity p^ takes 



values [ — J + 1 , 



•,L£J 



and N,j is the number of 



lattice points in the /1 direction. Thus, if we keep the 
physical size constant and indicate with p^ the mo- 
mentum components in the continuum, we find that the 
lattice components p^ = ap^ take values in the interval 
(— 7T, 7r] when a — > 0. 

In order to analyze the effects due to the breaking of 
rotational symmetry we consider two types of lattice mo- 
menta, i.e. 

• asymmetric: p\ — pi = p% = and pi =/= , 

• symmetric: p\ = p2 = P3 = Pi 7^ . 
For a symmetric lattice (L^ = L) we obtain 



P = Pi 
in the first case and 



2 sin 



P 



2p fl = 4 sin 



n pa 



7T pa 



(89) 
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in the second one. In both cases p takes the values —M- + 
1 , . . . , y since = N and N is always even. 

Let us note that the term contributing to |</> b (p)| when 
k = p is given by 



A. The ghost-gluon vertex 



Following the notation in Refs. 73, 90] , the ghost-gluon 
vertex function is given (on the lattice) by the relation 



6 50 Pp 4(0) G{p) 



(90) 



rr*(fe,p) 



V(A^(k)G cd (p)) 
[D(k))(G(S))(G(p): 



(94) 



Thus, if we neglect all the other terms appearing in Eq. 
(I57|) we find 



4>\p) | ) ex 



A b J0) G(p) 



[G(p)\ 



(91) 
(92) 



where s = k + p. ( D{k) ) is the lattice gluon propagator 
and A b (k) is defined in Eq. {73- At tree- level (on the 
lattice) one has [91| 



rn fc >p) = *gof bc %cos 



Thus, we can write 



(95) 



Therefore, if the ghost propagator has the infrared behav- 
ior (G(p)) ~ l/p 2+2K we obtain (in this approximation) 



T bcd {k, p) = i go f bc % cos ( ^ )T(k,p) (96) 



1 



3+2k 



(93) 



Clearly, for k = 0.5 one gets (G(p)) ~ 1/p 3 and 
( | </>(p) | ) ~ 1/p 4 - Let us recall that for the infrared ex- 
ponent k — in Landau gauge and in four dimensions 
- there are several analytic predictions (usin g D yson- 
Schwinger equations) Ha E3, E3, US 0, ES with 
values in the interval [0.52,1.00]. On the other hand, 
numerical studies have found a value of 0.25 — 0.35 
III HHH, both for the SU(2) and SU{3) cases. 
Recently, a value k = 0.22(5) has been obtained in the 
unquenched SU(3) case [52|. Note that these numeri- 
cal results for the infrared exponent n represent a lower 
bound, since the value of k increases when one is able to 
get data at smaller momenta using larger lattice volumes 
|57l l60l l87| . We investigate the infrared behavior of 4>(p) 
from our data in Section [V] below. There we also com- 
ment on the approximation obtained by neglecting terms 
with p =/= q in Eq. JS7J. 

As already stressed in the Introduction, the depen- 
dence of 4>(p) on A(q) [see Eq. l|63|l ] could explain the 
1/L 2 behavior observed for ( | <fi(p) \ ). On the other hand, 
as mentioned before, this result does not imply that <fi(p) 
is zero at infinite volume. 

Finally, from Eq. (|70|l we obtain that the off-diagonal 
symmetric components of (G b , c (p)) should also be non- 
zero if the off-diagonal anti-symmetric components are 
non- vanishing. However, in this case — due to the factor 
k^A^p — k)/p 2 — one should expect a more singular 
dependence on the inverse momentum 1/p and stronger 
finite-volume effects than those observed for the anti- 
symmetric components ( G b 7 c (p) ). 



and at the asymmetric point, with zero momentum for 
the gluon (k = 0), we find pi 



\)p 2 ^ 



/y„.V. ( .V - - l)p 2 ^ cos(^) 

Vf bcd (Al(Q) G cd (p)) 
(D(0))(G(p)) 2 ' 



(97) 



where we set E(p) = T(Q,p). Note that the combination 
f G (p)/N c is the same that appears in the definition 
of 4> b {p). This gives 



i 



90 (N 2 



l)p 2 



E 



cos 



/ 7rp M a 



V(A^0)^(p)\ 
(D(0))(G(p)y 



(98) 



On the other hand, since in the above equation the quan- 
tity (j) b (p) is multiplied by 4(0), one E e ^ s that the scalar 
function S(p) does not change sign under the transfor- 
mation (|73|l . This explains why there is no sign problem 
[73L l90| when evaluating the ghost-gluon-vertex renor- 
malization function Z^ 1 ^) — S(p). 

Let us also note that the factor cos (wp^ a/L^), which 
appears in the tree-level result l|95l) and in the denomi- 
nator of Eq. H98fl . is a discretization effect that goes to 1 
in the formal continuum limit a — > 0. Due to the relation 
between S(p) and (j> b (p) [see Eq. (|98|) ]. it is likely that 
similar discretization effects come up also in the present 
work. This will be verified in the next section. 
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FIG. 2: Histogram (solid curve) of <j> 3 {p) for the lattice volume 
8 4 with f3 = 2.4 and p = (0,0,0, 1), using 10,000 configura- 
tions. For comparison we also plot a Gaussian distribution 
(dot-dashed line) with the same mean value and the same 
standard deviation. 



V. RESULTS 

We start our analysis of the data by considering the 
two Fourier-transformed quantities 

Gb e d (p) = ^ E Qbd ^y) cos [p 0= - v) 1 (as) 

x y 

G h d (P) = ^E G M (x,y) sin [p (x -y)}. (100) 



As explained in Section IIIII we expect G bd {p) [re- 
spectively G bd (p)] to be symmetric (respectively anti- 
symmetric) in the color indices b and d. We have verified 
this by evaluating, for each configuration and for each 
momentum p, the quantities 



(101) 
(102) 



which should both be zero. We find that these two quan- 
tities are usually very small, being respectively of the 
order of 10 -3 and 10 -6 . At the same time, we checked 
the hypothesis of strong diagonal-dominance for the sym- 
metric components G bd (p). We obtain that the relation 




\G?(p)\ > 



E G 

d=£b 



bd 



(P) 



(103) 



is satisfied for almost all momenta p and configurations 
(except for 14 cases out of 280800). This is expected 
since, due to the global color symmetry and for a zero 
ghost condensate v, one should have (G e e d (p)) oc S ed , 
after averaging over the gauge-fixed Monte Carlo con- 
figurations. Thus, it seems reasonable to assume that 
G e e d (p) is strongly diagonally-dominant already for each 
configuration. This should also happen if the ghost con- 
densate v is nonzero and small. We also considered the 
propagator G bd {k 1 p) with k ^ p for the lattice volume 8 4 
and j3 = 2.4 with k given by (0, 0, 0, 1) in the asymmetric 
case and by (1,1,1,1) in the symmetric one. In this case 
we find that the relation 



\G?(k,p)\ > I 



E G ' 

d=tb 



hd 



(104) 



is satisfied in about 75% of the cases. 

We then evaluated the average of [see Eq. i|86fl] 



■ f bcd f bcd 

— G cd (p) = —G cd {p) 



(105) 



for b = 1,2,3. We find that this quantit y is c onsistent 
with zero, within one standard deviation |110| . Indeed, 
the statistical fluctuations are always very large, being 
larger than the central value for more than 60% of the 
data points and larger than 1/3 of the central value for 
all our data. On the contrary, when considering the ab- 
solute value of (j> b (p) we find a relatively small error, 




FIG. 3: Results for the Binder cumulant U [defined in Eq. |0] 
for the quantity (j> b (p) as a function of the lattice side L (in 
GeV - ) for all lattice volumes V, couplings f3 and momenta 
with p — N/4. We show the data corresponding to asymmet- 
ric momenta [for = 2.2 (□),£ = 2.3 (O) and f3 = 2.4 (A)] 
and to symmetric momenta (with the corresponding filled 
symbols for each (5). Errors have been estimated using the 
bootstrap method with 10,000 samples. 
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_ 10 
> 

(5 1 



0.1 
0.01 
0.001 
0.0001 
1e-05 



0.5 1 1.5 2 2.5 3 3.5 4 

p (GeV) 

FIG. 4: Results for 4> a (x), 4> m (□) and (j> r (o) [defined in Eqs. 
Hl()rifl - Hl()8|l ] as a function of the momentum p for lattice vol- 
ume V = 24 4 and f3 = 2.2. We show the data corresponding 
to asymmetric and to symmetric momenta. Note the loga- 
rithmic scale on the y axis. Errors have been estimated using 
the bootstrap method with 10,000 samples. 




0.5 1 1.5 2 2.5 3 3.5 4 

p (GeV) 

FIG. 5: Results for 4> a (x), <j) m (□) and <f> r (°) [defined in 
Eqs. H10ri^ - I|108[l ]. divided by cos (it p a / L) , as a function of 
the momentum p for lattice volume V = 24 4 and [3 = 2.2. We 
show the data corresponding to asymmetric and to symmetric 
momenta. Note the logarithmic scale on the y axis. Errors 
have been estimated using the bootstrap method with 10,000 
samples. 



of less than 10% of the central value in 99.7% of the 
cases. As discussed in the Introduction, this is a clear 
indication that the data for <fi b (p) do not correspond to 
a Gaussian distribution. We have verified this by com- 
paring a histogram of the data with a Gaussian curve. 
Then, as for the 3d 0(4) nonlinear a-model in the bro- 
ken phase (see the Introduction) , one finds (see Fig. [5J 
that the shape of the statistical distribution of the data 
is clearly not Gaussian. (However, notice that the shape 
of this non-Gaussian distribution is not the same as for 
the spin-model case.) We also find that working at con- 
stant physics, i.e. with approximately fixed lattice side L 
and momentum p (in physical units) , the deviation from 
a Gaussian distribution increases when the lattice cou- 
pling p increases. Finally, we have looked at the Binder 
cumulant U [see Eq. for the order parameter 4> b (p). 
We find that U is clearly different from zero for all lattice 
sizes N, couplings (3 and momenta p, and seems to con- 
verge to a value U ~ 0.45 as the lattice side L increases 
(see Fig. [3Jl . This limiting valu e is essentially indepen- 
dent of the momentum p |lllj . The data show good 
scaling for the three (3 values considered. As argued in 
the Introduction, the large fluctuations for <fi b (p) and the 
deviation from a Gaussian distribution could signal the 
spontaneous breaking of the global SL(2, R) symmetry, 
implying a nonzero value for the ghost condensate v. 

In order to plot the data for |</r(p)|, we average over 



the color index b considering the three quantities 

Mp) = ^E(I^)I) ( 106 ) 



= (JoEI^(p)P) (107) 



Mp) = W)i 2 >' ( 108 ) 

i.e. we take (respectively) the average of the absolute 
value of the components of </>, the average of the abso- 
lute value of <f> and the root mean square of <j). The 
results are shown in Fig. 0] In this figure the breaking 
of rotational symmetry is evident: the asymmetric and 
symmetric momenta clearly form two different curves. In 
particular, the three quantities decrease significantly at 
p sw 1.85 GeV in the asymmetric case and at p w 3.7 GeV 
in the symmetric one. 

The interpretation of the data is simpler (see Fig. 03) 
if one rescales <t> h (p) with the factor cos (up a/L) that 
appears in the lattice formula for the ghost-gluon-vertex 
renormalization function [see Eq. H98|l in Section TlV Aj . 
Moreover, as can be seen in Fig. the effects due to the 
breaking of rotational symmetry are relatively small. Fi- 
nally, we can multiply the data by L 2 to obtain the single 
curve shown in Fig. Again, we find good scaling for 
the three f3 values considered. Let us recall that, for the 
gluon and ghost propagators in the £{7(2) case \1% l88j . 
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> 10000 
CD 



1000 



O 
O 



100 




Data considered 


d.o.f. 


r 


2 


X 2 /d.o.f. 


All data 


208 


37.3(7) 


3.52(2) 


6.54 


V = 24 4 and /3 = 2.2 


20 


36(1) 


3.65(5) 


3.09 


All data 


96 


39.5(6) 


3.79(3) 


5.33 


and p < 2 GeV 










V = 24 4 , P = 2.2 


13 


36(1) 


3.80(5) 


2.50 


and p < 2 GeV 










All data, p < 2 GeV, 


24 


49(1) 


3.82(5) 


3.44 


only symm. momenta 










All data, p < 2 GeV, 


70 


37.7(4) 


3.80(2) 


3.35 


only asymm. momenta 










V = 24 4 , f3 = 2.2, 


9 


34.3(4) 


3.80(2) 


0.82 


only asymm. momenta 











p (GeV) 



TABLE III: The parameters r and z obtained by fitting 4> m 
to the formula 111211 . We considered different sets of data. 
The fits have been done using gnuplot. 



FIG. 6: Results for 4> m [defined in Eq. 11071 1. rescaled by 
L 2 / cos (irpa/L), as a function of the momentum p for all lat- 
tice volumes and P values considered. We show the data cor- 
responding to asymmetric and to symmetric momenta. Note 
the logarithmic scale on the x and y axes. Errors have been 
estimated using the bootstrap method with 10,000 samples. 



scaling has been observed in the range (3 S [2.15,2.8]. 
Note that, from our analysis in Section ITVl we would ex- 
pect (in the infrared limit) a behavior of the type G(p) /p 
for the three quantities defined above, if the contribution 
proportional to ^4(0) is dominant. 

Let us recall that in the SU(2) case the inverse 
Faddeev-Popov matrix — in the condensed Overhauser 
vacuum <fi b cx S b3 — should be given (in momentum 
space) by 0] 



G cd (p) 



p 2 S cd + v e cd 
p 4 + v 2 



(109) 



where e cd is the anti-symmetric tensor with c, d = 1, 2 
and v is the ghost condensate, given (at one loop and for 
N c = 2) by 







Xvac (2N f - 13) 


Po 




3 



With Xvac 

find 



(110) 



0.539 A|j~. Since in our case Nf = 0, we 



H 0.539 Afa 



2 34 A2 



(111) 



Then, with A^ » 250 MeV, we obtain v » 0.146 GeV 2 
A fit of the type 



cos 



( 7T p a \ 



r 
pz 



(112) 



works quite well for our data, giving for the parameters 
r and z the results reported in Table IlTTI Note that the 
exponent z increases systematically when one considers 
data at small momenta only, suggesting that z would 
approach 4 if one could really probe the infrared limit. 
We also tried a fit using the form r/(p z + v 2 ). The results 
obtained in this case are very similar to those reported 
in Table ITTT1 with v always consistent with zero within 
error bars. Furthermore, for the ghost propagator G(p) 
we find that a fit to r/p z yields z in the interval 2.33 — 
2.46, depending on the set of data considered. Thus, 
the quantity (f>(p) has a stronger infrared behavior than 
G(p) /p, which is obtained by neglecting terms with p =/= q 
in Eq. (JH3 . 

By comparing the above fitting formula to the theo- 
retical prediction (|109fl . i-e. by c onsidering (in the anti- 
symmetric case) the Ansatz |l!2j 



G cd ( P ) = 



rL~ 



we can state that 



v <§; r Li 



In our simulations we have 1/L max 
Pmin ~ 0.245 GeV. Thus, using r w 
we obtain our final estimate 

v < 0.058 GeV 2 , 



(113) 



(114) 

fa 0.039 GeV and 
38 (see Table ILTl. 



(115) 



Note that, since 38 ps (27r) 2 and p m in ~ 2ir/L max , the 
two inequalities in (|114(l give practically the same esti- 
mate for the ghost condensate v. One can also try a fit 
of the data using the function r/(p z + v 2 ) with a fixed 
small value for v . Then, for v < 0.025 GeV 2 the values 
of x 2 /d.o.f. reported in Table ITTT1 do not change signifi- 
cantly. 
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FIG. 7: Results for A a (0) (x), A m (□) and A r (0) (o) [denned 
in Eqs. pi6^ - l|118|l ] as a function of the inverse lattice side 
1/L for all lattice volumes and (3 values considered. Errors 
have been estimated using the bootstrap method with 10,000 
samples. We also report the fit of A m to the function b/L s 
with b = 10(7) and s — 2.4(2). When considering the data for 
A a {0), one finds for the exponent s the same value as above, 
while for A r (0) one gets s — 2.0(1). In all cases, the fits 
have been done using gnuplot and we consider only the data 
corresponding to 1/L < 0.07 GeV. 



Now, we consider the dependence of |A^(g)| on the 
lattice side L in physical units. In analogy with Eqs. 
(|106fl - H108fl . we define the three quantities 



A a (q) 



A m (q) 



A r (q) 



A h M i 




(116) 



(117) 



(118) 



As shown in Fig. [7\ we find that the (absolute value of 
the) Fourier-transformed gluon field at zero momentum 
A(0) goes to zero as L~ s with s rs 2 for the three quanti- 
ties de fined above. Similar results can be obtained when 
This is a confirmation that the finite-size 
effects observed for the quantity <fi(p) ar e indeed induced 
by A(q). We have also checked that the sign of 4> 3 (p) 
can be changed using the global gauge transformation 
reported in Eq. l(7^|) . i.e. by changing the sign of A^(q). 

Finally, we study the off-diagonal symmetric compo- 
nents G b , c (x,y) by considering the quantity rjr(p) — 
[G c e d (p) + G d e c (p)]/2, with b ^ c,d and b = 1, 2, 3. 
Let us stress that, even though these components are 




p (GeV) 

FIG. 8: Results for ip m [defined in Eq. H120H 1. rescaled by 
L 2 '° I cos (npa/L), as a function of the momentum p for all 
lattice volumes and (3 values considered. We show the data 
corresponding to asymmetric and to symmetric momenta. 
Note the logarithmic scale on the x and y axes. Errors have 
been estimated using the bootstrap method with 10,000 sam- 
ples. 



symmetric under the transformation (|73|l . they are not 
necessarily positive [see Eq. I|69fl]. Thus, also in this case 
we should consider their absolute value. In analogy with 
the analysis reported above for (f> b (p), we can evaluate 



\^{P) 



Mp) = (^Ei^)i 2 ) 
MP) = Jd E iw)i 2 ). 



(119) 



(120) 



(121) 



We obtain that these quantities produce a single curve 
when rescaled by L 2 - 5 / cos (irpa/L) (see Fig. [HJ- Note 
however that the rescaling does not work as well as for 
the ci's. Also, the rescaled data can be reasonably well 
fitted using the fitting function s/p z with z»5. This is 
in agreement with the observation made before Section 
IIV Al and in partial agreement with the results reported 
in Ref. H3- 

We have done the analysis reported above also for the 
gauge-fixed configurations obtained using the smearing 
method jT^. As for the ghost-gluon vertex jfj H3, 
we find that the contribution of Gribov-copy effects (if 
present) is always zero within error bars. In particular, 
the exponent z obtained when fitting the functions 4>{p) 
is close to 4 at small momenta also when using the smear- 
ing method. 
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VI. CONCLUSIONS 

We study numerically the off-diagonal components of 
the momentum-space ghost propagator G cd (p) for SU{2) 
lattice gauge theory in minimal Landau gauge. We see 
clear signs of spontaneous breaking of a global symmetry, 
using the quantity 4> b (p) — e bcd G cd (p)/2 as order param- 
eter. As in the case of continuous-spin models in the 
ordered phase, we find indication of spontaneous sym- 
metry breaking from two (related) observations: 1) by 
comparing the statistical fluctuations for the quantities 
4> b {p) and |</> b (p)|; 2) from the non-Gaussian shape of the 
statistical distribution of 4> b (p), which can be observed by 
considering a histogram of the data or by evaluating the 
so-called Binder cumulant. Since, in Landau gauge, the 
vacuum expectation value of the quantity <p b (p) should 
be proportional (in the so-called Overhauser channel) to 
the ghost condensate v [l3|. it seems reasonable to con- 
clude that the broken symmetry is the SL (2, R) symme- 
try, which is related to ghost condensation [y, |42] ■ Let us 
note that, from our data, the Binder cumulant U seems 
to be approximately null at small lattice volume and to 
converge to a value U rs 0.45 for physical lattice side 
L > 15 Gev -1 w 3 fm, corresponding to a mass scale of 
less than 100 MeV. 

As stressed in the Introduction, we have shown that 
the dependence of <j) b (p) on A d (q) [see Eq. {87)] can ex- 
plain the 1/L 2 behavior observed for ( | 0(p) | ). We have 
also shown that (j) b (p) has discretization effects similar to 
those obtained for the ghost-gluon vertex |jj |9(J • Then, 
using the rescaled quantity | L? 4> h {p)/ cos(irpa/L) |, we 
find a behavior p~ z with z « 4, in agreement with an- 
alytic predictions ^ is important to note that this 
result is not obtained just by considering the approxima- 
tion <j)(p) k, G{p)/p [see Eq. 11920 ]. since our data for the 
ghost propagator G{p) have an infrared behavior given by 
l/p 2A . We therefore view the similar momentum depen- 
dence obtained numerically and analytically for <fi b (p) as 



a further indication of ghost condensation. On the other 
hand, from our fits we find that the ghost condensate v 
is consistent with zero within error bars, i.e. the quantity 
| L 2 <fi b (p)/ cos(-Kpa/ L) | does not approach a finite limit 
at small momenta, at least for p > 0.245 GeV. Using 
the Ansatz (|113fl we obtain for the ghost condensate the 
upper bound v <C 0.058 GeV 2 . More precisely, our data 
rule out values of v greater than 0.058 GeV 2 w (240 
MeV) 2 but would still be consistent with a ghost con- 
densate v < 0.025 GeV 2 « (160 MeV) 2 . As said in the 
Introduction, in analytic studies 0, 0, |23|, |92 one finds 
that the ghost condensate induces a tachyonic gluon mass 
proportional to y^v, which modifies the dynamic gluon 
mass related to the gauge condensate (A^A^). Thus, 
one should expect a relatively small ghost condensate v 
in order to obtain a global (non-tachyonic) gluon mass. 
Let us recall that a dynamic gluon mass of the order of 
a few hundred MeV has been considered in several phc- 
nomenological studies [93, |9J, 19a . 19a. 1971] . A similar mass 
scale was also obtained in numerical studies of the gluon 
propagator in Landau gauge [H, |9^, llOOl Il0l| . 

We plan to extend our simulations to larger lattice vol- 
umes, up to 48 4 , which would allow us to detect a ghost 
condensate as small as v w 0.01 GeV 2 = (100 MeV) 2 . 
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